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SUMMARY 


The  objective  of  the  research  work  presented  in  this  re|ort  was  to  evaluate  existing 
finite  difference  methods  and  select  and  apply  the  one  most  appropriate  for  computing  the 
propagation  of  multiple  shock  waves  in  a  two-phase  rocket  combuston  chamber. 

In  order  to  be  acceptable  for  the  intended  application,  a  finite  difference  integration 
scheme  must:  preserve  the  high  frequency  content  of  the  waveforms,  be  relatively 
nondissipative  and  nondispersive  after  many  wave  cycles;  be  capable  of  describing  a  shock 
wave  as  a  sharp  discontinuity;  and  be  capable  of  properly  treating  the  reflection  of  shock 
waves  from  boundaries  and  the  partial  reflection  and  transmission  ot  discontinuities. 

1  he  finite  difference  schemes  of  the  Split  Coefficient  method,  the  X -scheme,  Rubin 
and  liurstein,  MacCormack,  Lax- Wendroff,  Godunov,  Rusanov,  the  Flux-Corrected-  I  ransport 
scheme  of  boris  and  tiook,  Chorin's  implementation  of  Climm's  Method,  the  Hybrid  scheme 
of  Harten  and  Zwas  and  the  Artificial -Compression  method  of  Marten  were  evaluated.  All 
of  these  schemes  are  finite  difference  approximat  ions  to  the  derivatives  arising  in  the 
conservation  laws  and  can  treat  an  arbitrary  system  of  conservation  laws. 

As  a  first  step  in  evaluating  the  potential  of  these  schemes,  they  were  utilized  in  the 
solution  of  the  shock  tube  problem  for  the  one-dimensional  tulerian  form  of  the  gas 
dynamic  conservation  equations  for  an  inviscid,  non  heat-conducting  fluid  and  for  the 
solution  of  the  linear  wave  equation  problem.  The  final  test  problem  was  the  solution  of 
the  one-dimensional  nonlinear  hyperbolic  equations  describing  finite  amplitude  wave  and 
shock  propagation  in  a  closed  end  tube.  In  addition,  the  ability  to  spectrally  analyze  the 
computed  results  was  developed.  This  capability  simplifies  the  interpretation  of  the 
complex  waveforms  and  facilitates  comparisons  among  the  various  finite  difference 
integration  schemes. 

The  results  of  applying  several  of  the  alxrve  mentioned  schemes  to  the  last  test  case 
are  presented  and  discussed  in  this  report.  It  had  Ireen  concluded  that  for  the  present 
problem,  a  method  based  upon  the  combination  of  the  I  ax -Wendroff ,  Hybrid  and  Artificial 
(.oppression  scheme  was  found  to  be  superior  to  the  other  schemes  tested.  I  his  sc  heme 
was  incorporated  into  the  nonlinear  instability  program  developed  by  I  evino  and  (  ulick  that 
described  the  comlxjstion  and  flow  inside  a  solid  rocket  nx>tor.  Finally,  the  ability  of  this 
scheme  to  treat  various  initial  disturbances  in  a  solid  rocket  motor  sucn  as  a  first  mode 
disturbance  and  standing  or  traveling  pulses,  is  demonstrated. 
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SIC  I  ION  1 


I  NT K unuc:  I  ION 

Ibis  report  presents  the  results  ot  cm  investigation  designed  to  select  a  satisfactory 
method  tor  computing  the  propagation  of  steep- f ranted  shocklike  waveforms  in  a  two-phase 
rocket  combustion  chamber. 

tactical  solid  rocket  motors  are  frequently  subject  to  a  combustion  instability 
problem  at  some  |X)int  in  the  design  cycle.  When  instability  is  encountereo,  it  can  take 
one  of  several  tor  ms,  e.g.,  linear  or  nonlinear,  longitudinal  or  tangential.  Over  the  last 
twenty  years,  considerable  resources  have  been  expended  to  understand,  predict,  control, 
and  eliminate  combustion  instability  in  solid  rocket  motors.  Most  of  this  effort  has  Ijeen 
devoted  to  linear  instability  problems,  and  as  a  result,  such  problems  can  now  be  treated 
in  a  rational,  cost  effective  manner.  In  comparison,  little  work  lias  been  accomplished 
towards  the  understanding  and  resolution  ot  nonlinear  combustion  instability  problems. 
Ihus,  when  nonlinear  instabilities  are  encountered,  the  solution  is  too  often  an  expensive 
cut -and -try  process. 

Linear  instability  is  characterized'  by  small  amplitude,  sinusoidal  oscillations  that 

originate  from  the  amplification  of  infinitesimal  random  disturbances  in  the  motor  chamlier. 

Nonlinear  instability  is  usually  characterized  by  large  amplitude  oscillations  having 

steep- f  routed,  shocklike  waveforms,  and  is  initiated  by  random  Unite  amplitude  events  such 

as  the  expulsion  of  an  igniter  or  insulation  fragment  through  the  nozzle.  Nonlinear 

instabilities  art*  modeled  using  lx>th  "exact"  and  "approximate"  mathematical  techniques. 

I  ne  "exact"  methods  of  Levine  and  Oulick,^  and  those  of  booker  and  / inn seek  to  solve 

inrrie  rical  ly  the  nonlinear  partial  ditlerential  equ.it  ions  governing  Ixrth  the  mean  and  time 

dependent  tlow  in  the  combustion  chamber,  as  well  as  the  comlmstion  res|jonse  ot  the  solid 

\  4 

propellant.  I  lie  "approximate"  methods  ot  <  ulick,  and  ot  hnvell,  et.al.,  utilize 


I.  I  evine,  ).  N,  and  i  ulick,  I.  I.  (  .,  ’Nonlinear  Analysis  ot  Solid  k<><  ket  (  ombust  ion 
Instability,"  AIKII  technical  Ke|>urt  IK-74-4S,  October  1*1/4. 

1.  booker,  I ).  I.  ami  /inn,  It.  I.,  "Numerical  Solution  ot  \xial  Instabilities  in  solid 
t’ropellant  Kocket  Motors,"  tilth  |ANN\I  Combustion  Meeting,  Yol.  1,  Naval  war  <  ollege, 
Newjort,  K.I.,  August  1 0 7  S ,  (  I’l  A  I’ub.  Z4 

t.  bulick,  t.  I.  "Nonlinear  behavior  < > t  Acoustic  Waves  in  < omt.ust n >u  i  hambers,"  totb 

)  \NNAf  t  ombustion  Meeting,  Vol.  I,  Naval  war  <  ollege,  Newport,  K.I.,  August  1**7  t. 

4.  (dwell,  i.  A.,  i'adi'iauaooan,  M.  s.,  and  /inn,  b.  I.,  ’  \pproximut->  Nonlinear  Vmihsis  ot 

Solid  Koc  ket  Motors  and  I  -  I'm  rue  r  s, "  AIKII  -  I  K  -  77  -  4<t,  |uh  1*177. 
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perturbation  techniques  anil  harmonic  analysis  to  reduce  tlie  "i  ivcrmuj;  <ii  t  lerent  tal 
equations.  I  ...  h  ot  tnese  methods  lias  certain  advantages,  disadvant  ages,  and  limitations. 
I  he  two  classes  of  approaches  complement  each  other,  and  etlorts  to  develop  troth  further 


are  warranted. 

I  lie  existing  "exact"  nonlinear  instability  programs  were  developeil  a!  tout  seven  years 
ago  and  are  not  capable  of  treating  the  multiple  shock  steep- 1 routed  ty|re  ot  iris!  anil  it  res 
that  occur  in  reduced  and  minimum  smoke  tactical  motors  developed  since  then.  I  he 
oo/ective  ot  the  present  research  is  to  extend  and  improve  the  "exact"  model  developed  by 
Levine,  et.  al.,  in  reference  1  to  the  [xhnt  where  it  can  be  used  In  motor  designers  as  a 
tool  to  aid  in  tlie  efficient  resolution  ot  such  nonlinear  longitudinal  instability  proble  ns. 
I  he  tirst  phase  ot  this  research,  the  results  ot  winch  are  re|x>rted  herein,  was  devoted  t<i 
improving  the  finite  difference  numerical  technique  used  to  solve  tlie  equations  governing 
wave  propagation  in  the  combustion  chamber. 

In  order  to  tie  acceptable  tor  the  intended  application,  a  unite  nit  lerent  <■  technique 
must  (a)  preserve  the  Itigli-f  requettcy  content  ot  tne  waveforms;  {b,  tie  relatively 
nondissipat  ive  and  non -dispersive  after  many  wave  cycles,  u,  be  capable  ot  bestriding  a 
shock  wave  as  a  sharp  discontinuity  without  generating  overshoots  or  undershoots  u|ion 
crossing  the  discontinuity;  and  Id)  ne  capable  ot  partial  reflect  ion  and  transmission  at 
area  discontinuities.  It  should  be  |mint<  d  out  that  in  solving  a  combustion  instability 

problem,  numerically  induced  pro-  and  post-shock  "wiggles"  an  not  gist  impair  the 

accuracy  ot  the  solution  but  could  also  falsely  "trigger"  nonlinear  instabilities  anu  tone 

the  scheme  to  pick  a  nonphysical  solution. 

It  is  very  difficult  tor  any  single  finite  difference  scheme  to  satisty  .ill  ot  the 

aforementioned  requirements  simultaneously .  lor  example*,  several  artificial  viscosity 
schemes  have  been  developed  to  damp  pre-  and  |iost- shook  oscillations.  ■  inwove t,  su<  ii 
artificially  introduced  (.fit  fusion  also  smears  out  the  discontinuity  and  eventually  damps  the 
high-frequency  modes  that  are  part  ot  the  physical  I  low  Mold  inside  the  motor.  Moreover, 
such  artificial  damping  can  lie  comparable  in  magnitude  to  the  usual  net  gams  or  losses  ot 
acoustic  energy  111  roi  ket  motors,  nonce.  Its  presence  would  seriously  limit  the  validity  ot 
the  results.  In  addition,  use  ot  an  artificial  viscosity  would  hamper  eltorts  to  determine 

>.  I  ovine  and  <  ill  i<  k,  <  ;.  t  il . 

li.  I  apldus,  \.,  "  \  Molar  lied  Mum  k  I  air  illation  by  Second  1  'rder  I  mite  l  •  1 1  |e  rein  es ,"  |.  ot 

y  omp.  Iliys.,  \nl.  2,  pl>.  lr.4-17,-,  1  ‘ill/. 
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tiu-  actual  damping  ot  tiits  pli.is  -  11  '...ms  in  metallized  solid  pro|>ellaiit  rocket  motors. 

In  this  am"  1  :  ...u,  it  should  or  (Minted  out  that  shoe  k  -  t  itt  mg  schemes  that  t  rent  1 1 it ■ 

shock  as  an  internal  Uiuntiary  are  impraet  ical  tor  tins  appl  it  at  ion,  tine  to  the  large  iui.ia.er 

ot  shocks  traveling,  interacting,  anti  reflecting  uisitie  varialile  t  roSs -set  t  ional  area  motors. 

'sum I j r I v  ,  t  mitt*  tiitterence  schemes  (such  as  the  '-scheme  or  the  spl it  - 1 net  t it  iei:t 

a 

.Method  )  that  are  modeled  alter  anti  exploit  the  mathematical  theory  ot  the  methoii  ot 
characterist  ics  are  impractical  tor  our  specitic  problem.  shock -capturing  imp  licit 

tiitterence  schemes  otter  no  particular  advantage  tor  the  present  prohle  e,  since  the 
physical  problem  ot  interest  requires  time  resolution  consistent  with  the  stability 
restrictions  ot  explicit  methods. 


Six  or  seven  years  ago,  the  task  ot  1  Hiding  a  suitable  finite  dittereme  scheme  .voulii 

have  lieen  virtually  impossitile,  as  then  existing  methods  such  us 
10  ,  I  !  Id 

'  I  act  orinack ,  I  ax  -  weilurot  t ,  cardunov,  etc.,  lacked  the  required  qualities.  Suite 

then.  In  >vx  ever,  several  special  schemes  designed  to  achieve  si  lock  resolution  without 

‘wiggles’  were  developed.  Among  them  are  the  I  lux-t  orrecteo-  I  ransjxtrt  scheme  ot  «.ris 
1  '  14 

and  book,  l  horin’s  implementation  .it  (jlimm's  method,  the  i  pst  ream  (entered  scheme 


7.  Moretti,  l,.,  "  I  lie  titmice  ot  a  line  I  V|H>ncient  loc  liniquc*  in  t  ,as  !  k  namies,"  I’lli-M 
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springe  r -A  e  rlag ,  New  AuTh,  I‘i71. 
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of  Van-Leer,'^  t fie  Hybrid  Scheme  of  Marten  and  /was,"'*’  and  the  Artificial  Compression 
1  7 

Scheme  of  Marten.  All  of  these  schemes  are  finite  difference  approximat ions  to  the 
derivatives  arising  in  the  conservation  laws  and  can  treat  an  arbitrary  system  oi 
conservation  laws.  In  addition,  they  are  all  fixed  arid  methods  arid  can  automat  ically 
handle  interactions  between  waves  of  different  families.  hi  order  to  evaluate  the 
suitability  of  the  aforementioned  techniques  for  the  present  purjx>se,  their  ability  to  treat 
a  number  of  simpler,  but  related,  problems  was  examined.  I  he  results  of  this  investigation 
are  summarized  in  .Section  3.  I  he  finite  difference  scheme  judged  most  promising  was 
then  incorporated  into  the  nonlinear  instability  program  described  by  Levine  and  Culick^1. 
The  results  of  a  number  ot  nonlinear  instability  solutions  are  presented  to  demonstrate  ttie 
effectiveness  ot  the  new  technique. 


15.  Van  Leer,  If.,  ’lowards  the  Ultimate  Conservative  Difference  Scheme  III, 
Upstream-Centered  Unite-Hit  ference  Scheme  for  Ideal  Compressible  Hows,"  j.  ot  Comp. 
Physics,  Vol.  3,  1977,  pp.  2(>  5-275. 

Id.  Marten,  A.  and  /was,  C.,  "Self  Adjusting  Hybrid  Schemes  for  Shock  Computations,'  |. 
of  Comp.  Physics,  Vol.  9,  197  2,  pp.  5P8-5H3. 

17.  Marten,  A.,  “  I  he  Artificial  Compression  Method  for  Computation  of  Shocks  and  (ion  tact 
Discontinuities:  III,  Self  Adjusting  Hybrid  Schemes,'  AiOSK  technical  Kojxirt  I K  -  77 -l)(>59, 

March  1977. 


1 1I.  Levine  and  Culick,  Op.  Cit. 


MCI  ION  2 
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2.1  sihk:k  iuijl  pkobllm 

As  a  first  step  in  evaluating  the  various  finite  ilifterence  schemes,  they  were  used  to 
solve  the  shock  tube  problem  for  the  one -dimensional  Lulerian  form  of  the  gas  dynamic 
conservation  equations  for  an  mviscid,  non -heat  -conifuct  ing  fluid.  ihese  solutions  were 
used  to  rate  the  difference  schemes  based  upon  such  criteria  as  resolution  of  the  shock 
(i.e.,  number  of  mesh  points  needed  to  describe  the  shock  discontinuity),  diffusion  and 
smearing  ol  the  shock  and  the  contact  discontinuities  with  time,  stability,  effect  of 
Courant  number,  and  computation  time. 

1  he  results  of  the  shock  tube  test  case  demonstrated  the  superiority  of  the  recently 
1  q-2  f  24  2r> 

developed  techniques  over  earlier  methods.  '  based  on  the  previously  mentioned 

criteria,  the  best  method  was  a  combination  of  three  technique's.  I  fie  basic  sclieme  of 

2 1>  2  7 

Lax-Wendroft  was  combined  with  a  llydrid  scheme  and  the  Artificial  Compression 

2  a 

Method.  this  method  (LW+I  f+ACM )  was  capable  of  producing  very  sharp  shock  waves, 
without  pre-  or  (lost -shock  wiggles. 

2.2  I.INLAK  WAVl  fcQUAHON 

I  he  various  finite  difference  methods  were  also  employed  to  solve  a  second  problem, 
i.e.,  the  linear  wave  equation.  The  results  were  used  to  evaluate  the  relative  diffusive  anu 
dispersive  errors  of  the  schemes  tor  harmonic  standing  wave  propagation  after  many  wave 


lb.  Boris,  Op.  (lit. 

20.  Chorin,  Op.  Cit. 

21.  Van  l  eer,  Up.  Cit. 

22.  ilarten  and  2 wa s ,  Op.  Cit. 
2  5.  Marten,  Op.  Cit. 

24.  MacCormack,  Op.  Cit. 

25.  Oodunov,  Op.  Cit. 

2t>.  lax  and  Wendroff,  Op.  Cit. 
27.  Marten  and  /was.  Op.  Cit. 
215.  Marten,  Up.  <  it. 
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cycles.  Ihis  problem  was  also  used  to  assess  tbe  ease  ot  imple.nentat ion  ot  lioundary 
conditions  tor  the  techniques  considered. 

I  lie  linear  wave  equation  solutions  demons!  rated  that  the  magnitude  ot  the  dittusive 
and  dispersive  errors  varied  significantly  from  one  method  to  another.  I  lie  l»V*tl  +  .\vM 
combination  method  produced  excellent  solutions  to  this  problem,  with  little  error  evident 
after  many  wave  cycles. 

2.3.  WAV  I  t’KOI'ACAIION  IN  A  ClOStl)  IUIU 

I  he  final  test  leading  to  the  selection  ot  a  method  to  be  utilized  in  the  nonlinear 
instability  propram  svas  the  solution  ot  the  one-dimensional,  nonlinear  hvperlxdic  equations 
describing  finite  amplitude  wave  and  shock  propagation  in  a  closed  end  tube.  Ii\ 
establishing  different  initial  conditions,  it  was  possible  to  compare  results  obtained  with 
the  various  techniques  tor  problems  ranging  from  almost  linear  (small  initial  disturbances j 
to  highly  nonlinear  (large  disturbances)  conditions,  over  .1  large  lumber  of  wave  cycles. 
In  addition,  results  obtained  with  these  techniques  were  compared  tor  problems  such  as 
interactions  of  different  wave  families  and  the  effect  ot  change  in  grid  size  and  (  ourant 
lumber  upm  the  t  requency  content  of  the  scheme. 

Some  of  the  results  of  the  last  test  problem- -  finite  amplitude  wave  in  a  closed 
tube- -will  be  presented  to  illustrate  souk*  of  the  differences  between  the  techniques.  As 
a  result  ot  the  numerical  error  associated  with  finite  difference  methods,  each  technique 
acts  as  a  numerical  filter.  I  he  ‘filtering’  effect  of  each  of  the  techniques  is  a  different 
function  of  frequency,  mesh  size,  etc.  lo  enhance  the  ability  to  discern  these  differences, 
the  results  of  the  test  problems  were  spectrally  analyzed. 

figures  1,  2,  and  3  demonstrate  the  effect  of  mmcrical  technique  on  the  time 

(■volution  of  acoustic  pressure  amplitude  .it  an  end  ot  a  closed  tube.  I  he  solutions  were 

initiated  with  a  first  harmonic  standing  wave  perturbation  having  an  amplitude  ot  20  „  ot 

the  mean  pressure.  1  he  numerical  schemes  used  in  the  figures  were  the  Iw+ll+M  \t 

method  (fig.  1),  the  Mact.ormack  scheme  (fig.  2)  and  the  Rubin  and  Hurstem 
29 

scheme  (lig.  5).  the  evolution  of  |xist-shock  wiggles  into  erroneous  higher  harmonii  s  is 
evident  in  tbe  latter  two  cases.  Spectral  analyses  of  the  results  shown  m  I  igures  1,  2, 
and  i  are  presented  in  figures  4,  r>,  and  l>  respectively. 


29.  Rubin,  !.  I.  and  Hurstein,  s.  /.,  "Difference  4ct  hods  tor  the  Invisiid  and  Visions 
tquations  of  a  C  ompressible  (.as,"  (ournal  ot  <  oipMit.it  tonal  I’liysn  s,  Yn|.  2,  19t>7,  pp. 

17  «-•»(,. 


igure  2  a).  1  me  I  volution  of  Pressure  Oscilllations  at  an  r.nd  of  a  lube  (MacCormack) 


H|jme  2(t>)  and  (c). 


IxpamW  Views  <d  the  1  lino  l  volution  of  the  I’ressun 
0s<  illations  at  an  t  nd  ot  a  I  ul«'  ;Ma<  (  ormai  k  ). 


T  j  .  I  rm-  Isolation  of  ’’’res'-ur**  Oscillations  .it  an  l.nd  of  a  Tube  {Rubin  and  Hurstein) 


figure  i(b)  and  (c). 


fcxpanded  Views  of  the  lime  I. volution  of  Pressure 
Oscillations  at  an  I  nd  of  a  I ube  (Rubin  and  Hurstein) 
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V unction  of  frequency  (Kubin  and  "urstein) 
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V  ijjure  b(e). 


I  irne  I  volution  of  Accumulative  I’M)  as  a  I  unction  of  Mode  Numlier 
(Rubin  and  Hurstem). 


With  tin'  I.A  +  i  I  +  \l  \l  technique,  the  initial  ■>  int*  wave  ill  .turbance  develops  into  ,i 
sharp,  triangular,  shock -typo  wave  torm  ( I  ig.  1;.  Ihis  is  the  expected  result  based  on 
the  closed  torm  analytical  solution  to  this  problem,  according  to  Morse  and  I  ( i  rcJ. 
f  igures  4a  to  4d  present  the  tunc  variation  ot  the  power  spectral  density  contained  m 
each  mode.  Ihese  figures  indicate  that  higher  inodes  up  to  the  1  1th  are  excited.  figure 
4e  shows  the  time  variation  ot  the  accumulative  power  spectral  density  and  indicates  that 
tor  this  numerical  scheme  there  are  no  erroneous  shitts  ot  power  spectral  density  in  (lie 
high  order  mixles.  I  he  acoustic  energy  distribution  varies  little  with  time  .once  a  shock  is 
formed;,  as  expected.  I  tie  transfer  ot  energy  from  the  fundamental  into  the  higher  mixles 
results  from  the  physical  process  ot  wave  steepening.  there  is  no  physical  process  tor 
transferring  energy  the  other  way,  i.e,  from  higher  inodes  into  the  iunoainental. 

figure  2  shows  the  time  evolution  <>t  acoustic  pressure  amplitude  at  an  end  ot  a 
closed  tube  when  utilizing  MacC.ormack's  scheme.  \n  expanded  view  ot  the  acoust  ic 
pressure  amplitude  between  nondimensional  times  ot  zero  and  fifteen  is  shown  in  figure  Jb. 
I  he  development  of  the  initial  ;x>st  shock  "wiggle,"  that  appears  after  the  first  cycle,  into 
erroneous  higher  modes  is  evident.  An  expanded  view  ot  the  acoustic  pressure  amplitude 


net  ween  nt  >nd  irnen  s  iona  1 

times 

5(J  and  (>D 

( f  ig. 

2c)  shows  that  the  ac 

oust  ic 

pressure*  wave 

form  that  at  earlier 

t  lines 

contained 

11-12 

harmonics,  has  now 

been 

reduced  to  a 

waveform  composed  ot 

4-5 

harmonics 

only. 

Ihe  absence  of  tile 

higher 

harmonics  is 

evidenced  by  the  discrete  humps  in  the  waveform.  litis  trend  is  also  substantiated  by  the 
spectral  analysis  of  tins  solution  (shown  in  figs.  r>a  through  ">e).  I  igure  lib  indicates  an 
excessively  high  percentage  of  acoustic  energy  in  the  eighth  and  ninth  harmonics.  At  the 
end  ot  the  test  (after  .SO  wave  cycles),  it  is  shown  ( I  ig.  a d j  that  tills  excessively  high 
percentage  ot  acoustic  energy  has  been  transferred  to  the  littli  and  sixth  harmonics. 
Since,  as  previously  mentioned,  there  are  no  known  physical  processes  that  can  cause  such 
a  transition  of  energy  from  a  higher  mode  to  a  lower  mode,  this  phenomenon  is  apparently 
a  result  ot  numerical  error.  figure  so  is  another  way  of  demonstrating  the  same  spurious 
result.  Hence,  the  time  variation  in  the  accumulative  power  spectrum  is  a  result  ot  ttie 
dispersive  error  ot  the  numerical  scheme  which  causes  pressure  signals  to  travel  at  the 
wrong  speed. 

I  he  results  obtained  utilizing  the  isiibin  and  liurstein  scheme  (figs  f  and  »>!  are 
generally  similar  to  the  results  obtained  by  utilizing  Mact  ormack's  scheme  i  I  igs  1  and 

td.  Uor>e,  I’.  M.  and  Ingard,  h.  V.,  Iheoretical  Ac  iiistus,  McCmaw  llill  book  (om>ativ, 
New  York,  Ibbli. 

t 


\  compari-ain  ot  the  results  indicates  that  tin*  tirst  jxist -shoc  k  wiggle  appears  after  the 
third  wave  cycle  (K.b)  compared  to  the  tirst  wave  cycle  (with  War  (  orma*  k ).  Also,  the 
percentage  ot  energy  contained  in  the-  t  undamental  mode  is  greater,  and  the  erroneously 
hir;h  energy  in  the  higher  modes  is  somewhat  lc*ss  with  the  Kudin  and  liurstem  scheme.  In 
tills  connection,  it  should  lie  mentioned  that  very  similar  results  to  those  obtained  by 
utili/iily,  tile  Kudin  and  liurstem  scheme  were  obtained  with  the  I  ax  -  Wendrof  t  scheme. 

I  he  accumulative  percentage  of  [lower  spectral  density  contained  in  the  respective 
harmonics  tor  the*  nondimensional  time  ot  7.2  to  14. ltd  is  shown  in  I  i^ure  7  tor  live*  ot 
the  numerical  schemes  considered.  I  he*  superiority  ot  the  I.V+II  combination  scheme  over 

the  standard  schemes  is  evident.  I  In*  addition  ot  the*  artificial  compression  to  this 

combination  serves  to  sharpen  shock  transitions  by  exciting  higher  harmonics  that  were  over 
damped  by  the  addition  ot  the  Hybrid  scheme.  Ilius,  compared  to  the  I  .V  and  Hyorid 
result,  the  l.V+lt  +  .UW  solution  has  slightly  less  energy  in  the*  fundamental  mode  and 
slightly  more  energy  in  the  higher  harmonics. 

All  of  the  earlier  methods,  i.e.,  WacGirmack,  I  ax-VVendrott,  and  Kuliin  and  liurstem 
were*  used  without  adding  an  artificial  viscosity.  1  lie*  effect  ot  an  artificial  viscosity  on 
the  solution  was  investigated  usinn  Hyman’s''  prc*riictor -corrector  scheme.  Kesults  utilizing 
this  method  with  two  different  artificial  viscosity  ciH'tl  ic  ie/)f  s  are  presented  in  I  inures  II 
and  d.  I  he  addition  of  an  artificial  viscosity  to  a  numerical  scheme  was  conceived  .is  a 
way  to  damp  post-shock  oscillations.  Artificial  viscosity  reduces  |x>st  shock  oscillations  at 
the  expense  ot  tiie  higher  harmonic  component s  of  the  waveform.  I  ipure  H  shows  what 
happens  when  a  high  value  ot  artificial  viscosity  is  employed  (A  equal  unity  in  Hyman’s 
method).  In  this  case  the  hi^h  artificial  viscosity  prevents  a  shock  from  ever  tormina  and 

the  deviations  from  a  perfect  sin**  wave  are  never  lar^e.  I  he  results  ot  the  spectral 

analysis  ot  this  solution  (fins,  till  anti  He)  show  the  absence  ot  higher  harmonic  content. 
Keduc  iny>  the  artificial  viscosity  coefficient  (  <S  =0.1,  tin*  lowest  value  at  which  Hyman’s 
method  remains  stable)  yields  a  much  steeper  waveform  (li^.  ha),  hut  one  whose  higher 
harmonic  content  is  still  less  than  it  should  lx*  (figs.  Oh  and  dc).  As  time  passes,  the 
action  of  tht*  artificial  viscosity  continues  to  damp  preferent tally  tin*  higher  harmonics, 
causing  the  solution  to  further  degenerate. 


11.  Hyman,  lames  M.,  "On  Kolxist  and  A*  <  urate  ‘.tetlxxls  tor  the  (  a  I*  ul.it  ion  ot 
Compressible  I  luiil  I  lows,’  t’art  I  to  be  published. 
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As  mentioned  earlier,  artificial  viscosity  schemes  would  hamper  efforts  to  determine 
the  actual  damping  of  gas  phase  oscillations  in  metali/ed  solid  propellants.  lo  illustrate 
the  similarity  between  the  effects  of  particles  and  artificial  viscosity,  the  previously 
described  closed  end  tube  problem  was  modified  by  the  addition  of  2'.,,  <  particle  to  gas 

weight  flow  latio)  of  5  micron  particles.  I  he  computed  results,  obtained  with  the  Kubin 
and  iiurstein  scheme,  are  shown  in  figure  10.  Comparing  figure  10  to  figure  1  shows  that 
the  initial  |)ost -shock  wiggles  (tig.  lOh  is  an  expanded  view)  did  not  develop  into  discrete 
humps  in  the  waveform.  As  expected,  the  spectral  analysis  tor  the  particle  case  did  not 
show  the  high  erroneous  spectral  content  in  the  higher  harmonics.  With  particles  present, 
the  damping  of  the  high-frequency  content  of  the  waveform  is  due  to  a  real  physical 

process.  this  is  in  contrast  to  the  similar,  but  nonphysical,  action  ot  an  artificial 

viscosity. 


f  me  [volution  of  the  Normalized  Oscillatorv  Pressure  at  >n  End  of  i  '  dosed  I  ut»e  filled 
with  Particles  tKiibin  and  Iturstein  Scheme). 
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Figure  10(h)  and  (c). 


txpanded  Views  of  the  Normalized  Oscillatory  Pressure 
Amplitude  at  an  f nd  of  the  C  losed  lube  filled  with  particles 
(Rubin  and  Iturstein  Scheme). 
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MOIOIx  SUUIIIIK". 

Iiu“  test  problem  solutions  we  re  used  to  rank  tin-  various  methods  aiioroutj;  to  the 

previously  mentioned  criteria.  It  was  com  luded  tli.it  tin*  l  ax  -  Wendrot  l  *  ■  r\  t  >  r  i<  i 

\ r 1 1 1  If  it* I  <.  (impression  scheme  yielded  tin1  l«‘st  results,  thus,  this  si  heme  v.as  inr  o r  |  m r at  eu 
into  the  nonlinear  instability  pro); ran;  uesi  r  iheil  l>\  I  evine  aim  t  ulnk. 

I  he  i|uas  i  -  one  -  dlinens  lona  I ,  Kvo -phase,  ei|uations  ot  i  ml  toll  t  >  lie  solve.,  nr.’  shown  in 
Injure  II  ill  lonservative  tone.  I  he  propellant  properties  anil  motor  ,, comet  r\ 

i  v  I  lllitr  u  al  ly  pertorateii  strain  ‘ih./i  ill  l.s  unties  loiip  >  were  taken  Iroo  I  evine  anil 
(  ul  it  k  *  '  to  taiilitate  i  oinp.t  r  isi  hi  svitn  the  earlier  results,  i.elore  presenting  the  results  .  1 1 
t  lie  motor  solutions,  it  shoiiln  tie  mentioned  that,  tor  tin-  motor,  |  to|H-llanl  a, in  ■  .ik' rat  in,, 
londitioiis  utili/ed,  a  linear  st.ilhlitv  analysis  shows  the  tunoamental  uimie  to  ih-  unstanle, 
while  all  ot  the  lu^ner  moiies  are  stable.  It  should  also  pe  |«unted  out  tiiat  at  the 

present  tune  velocity  1 1  >i  i[  >1  il>K  etteits  are  not  liu  limed  in  the  nonlinear  instability  . . . 

t  h.ure  1  Ja  sliovvs  the  pressure  tune  history  at  the  head  enu  ot  ttie  motor  i.ih  nlalm: 

us  in),  the  l,y+ll+U\t  scheme  and  an  inuial  tirst  mode  cisturh.inie  ampl  ituoe  ot  o.-l  ot  the 

mean  pressure.  Ihe  early  history  trom  t  -  d  to  1  >  nonoiiiension.il  time  is  shown  n 

expanded  torm  in  I  i);ute  12li.  Ihe  initially  sinusoidal  wave  is  seen  !■>  almost  imuiedi.il e h 
under, ;o  transition  to  a  sharp,  osi  illation  tree,  shnik-type  wave  to  rim  Iipure  I  .’i  st.ow  . 
the  pressure  histories  imth  at  the  head  end  anil  renter  ot  hue  motor  tor  the  time  niter,  r 
t  III  to  till.  In  this  time  interval,  the  wave  le.it  lies  a  limiting  amplitude.  1  y. imin.it  i,  ,n 

ol  the  phase  ditlereiu.  s  between  the  pressiiii'  at  oillerent  lor  itioiis  iiihh  ,ites  that  trie 

solution  is  neither  a  standm;;  wave,  not  a  Iravelliij;  wave.  Put  i  to  ilu  i, it  ion  ot  with. 

Ihe  same  profile  n  was  solved  usrip  the  i.ulnu  arid  hurstein  si  m  i  ie  fie  '<  t  io>! 
employei!  in  the  original  program  hy  I  evine  and  (  uln  k.  1  tallies  I  ta  s  and  ,  sn.isv  the 

head  end  pressure  history  and  illustrate  the  develo|«uelit  ot  the  initial  ,kis!  -  sh.  ,.  I  ,\  i  ;,de 

into  .ill  erroneous  setoiid  mode  osi  illation. 


t  J.  I  evine  and  t  ulnk,  '  ip.  <  it . 
it.  1 1  nd. 
hi.  hud. 
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NORMALIZED  PRESSURE  PERTURBATION 
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1 2  !'  )•  I  Views  of  the  lime  tvrilutiou  of  the  Norm.tli /<■<) 

Dsi  ill.it, >rv  I’rcsMin  \mplitude  Ht  the  lie. id  I  nil  «>t  the  M 
1st  IhirmoiiK,  St.indmn  IV.ive,  '  I’  -  0.4  : 1  v\  ♦  1 1  ♦  \(  M 
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N0RMM.I7.KD  PRESSURE  PERTURBATION 


1  he  marked  improvement  obtained  with  tin*  n<  w  solution  ><■  heme  lor  t  i  1 1 


V  \  1  1 . 


case  a. 1  x  ,4  r.it  1 1  v  1114. 


lu  examine  the  etlect  ot  initial  disturbance  amp  I  it  ude,  the  same  prune- n  a, is  mii.'-i. 
uxina  ttu*  I  >\  + 1 1  ♦  V. v,  v  heme  ami  initial  t  ( r  •>  t  mode  disturbance  •!<  1  ip  I  itutles  ot  •  amt  ...<■ 


ot  t  lie  me 

an  prexxule. 

!  tie  lie 

ad  end  pressure 

hist-cries  tier  these  c 

uses  .ire  n!m»»\ii 

1  14U  re  s  14 

and  1  >  r.-xj 

n-c  t  is  els  . 

In  all  three  1 

axes,  <1.4,  I) .J,  and  1 ) .< 

i  1  nut  i«i !  rnpl  it  mi 

shown  ri 

1  1411  rex  |  J , 

14,  and 

1  1  respectively. 

the  xo| ut  |. in  reac  hex 

ttu*  s.i’in*  !  rint  m 

amplitude  , 

it  al>  nit  1.  1  J 

ot  the 

m.-an  pressure. 

Ihix  result  was  quite 

■  surprising  in  vh*. 

1 )  t  thi1  ri*sults  iis  I  evinc  and  t  thick  which  demoiist  rated  <1  strong  ettett  it  initial 
disturbance  on  limit  my,  amplitude.  Hie  explanation  appear-,  to  lie  in  tne  tat  t  that  the 
present  solutions  .sere  tor  ^a-  -1111.4  while  the  earlier  solutions  At-re  tor  the  same  -i-itoi, 
operating  conditions,  anil  pro|iellant ,  but  with  I  ">  ,  ahminum  atkietl  ,d.ln  particle  ^as  weight 
t  lov\  ratio).  Hence,  it  seems  that  nonlinear  particle  damping  ettects  plav  a  eritic.il  role 
in  establishing  limiting  ampl  itude. 

I  he  spectral  analysis  results  cor  res|x>ndiii);  to  the  il.it1>  amplitude  case  are  shown  n 
I  14u.es  I  t>a  through  Inti.  these  results  clearly  demonstrate  the  excitation  ot  higher 
harmonics  -is  the  origin, il  dis turlfuncr  /;r o.vs  and  ste<*pens. 

■\n  important  application  ot  nonlinear  instability-  analysis  1-  the  prediction  ot  motor 
response  to  pulse  type  disturbances.  I  lit-  I  evine  and  t  iilick  ’  nonlinear  instability  pio^ra  n 
has  ,1  li  nitetl  abil-tv  to  generate  pulse  type  xvavtMorms  -  it  will  he  e\pand<*<i  soon  .  Ihis 

capability  was  used  to  test  the  ability  ot  the-  I  .'vi-ll+At  \l  method  to  trc>at  such  problem-. 
I  he  r  •  *  su  1 1  s  ot  t  .so  solutions  .sill  he  discussed.  In  noth  rases,  the  initial  pressure 

disturbance  .vavetorm  was  taken  to  be  ot  the  torm  sin  1  \.-i)  with  an  amplitude  ei|ual  to 

U.4  ot  the  mean  pressure.  I  his  produces  the  centered,  symmetric  waveform  shown  in 
I  1411m  I  ~ .  I  lie  ditlerence  nelAn-n  the  two  cases  was  the  initial  velocity  at  t  -  (I.  In 
one  i  use,  the  velor  ity  was  taken  to  he  !’/  sshilc-  in  the-  second  <  ase,  the  velor  its 
(nondiniensiona!  ,,ax  taken  t-i  he  zero.  I  he  tirst  case  represents  a  traveling  pulse.  I  he 

value  ot  veloc  ity  I-.  iakc'ii  sue  It  that  the  initial  pulse  represents  a  risht  traveling  wave.* 

I  lie  second  1  .nr  (  or  responds  to  a  “  standing"  pulse.  I  he  pulse  propagates  as  t  hr  -  sum  id 

equal  left  and  rp;ht  t  ravelin);  waves,  each  bavin)*  halt  the  initial  amplitude. 

1  >.  I  Old. 

li,.  loni. 

•Actually  sett  in);  V  -  ’I’/  :  only-  produces  a  pure  and  rq.lit  tr.iveliii);  was*-  111  tne  linear 

I  unit  as  1  I’  approac  lies  I). 
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I  hi1  calculated  pressure  lilstor  es  at  the  head  end  ot  the  motoi  tor  eai  h  ot  these 
ilisturuaiices  are  shown  in  figure  "  ;i  y  traveling)  and  figure  In  stanuinK  •  I  in-  Uranutn 
difference  between  tin-  results  demonstrates  the  iinportaiice  ol  tlie  velo<  its  (iistutbame 
associated  with  <i  pressure  pulse.  I  he  traveling;  pulse  IS  Iininedi.it  el\  t  ruts?  rined  into 
st  ee(>  - 1  rout  <“d  shock  type  waveform  and  decays  until  it  reaches  the-  same  limiting,  implitade 
as  the  solutions  started  from  first  mode  sinusoidal  dist  uriijni  c>s.  Spectral  aiial.ses  .it  this 
solution  tor  three  time  intervals  are  she iwi  in  I  inures  ilia  to  JUi  .  \t  e.irlv  times  lip. 

20a)  the  traveling  pulse  disturbance  contained  i  lur^e  percent  ,i);e  cd  the  tundament.il,  ml  a 

significant  higher  harmotic  content  is  also  evident.  \s  the*  limiting  amplitude  was 
approached,  at  later  times,  thc>  percentage  <*t  energy  in  the  tundamental  me  tisisccl. 

I  he  |>ressure  history  ot  the-  standing  pulse  disturbance  is  shown  in  I  incite  I  'i.  i  he 
time  variation  c  t  the  waveform  is  <|uite  complex  in  tins  rase.  I  he  spectral  mahsis 

results  shown  in  I  inures  J  la  to  2  1c  help  to  clarity  what  is  happening.  \t  early  times,  a 

symmetric  standing  pulse  centered  in  the  motor  contains  essentially  only  even  harmonic  s 
( I  ig.  2  1a),  with  the  second  harmonic  dominating;.  I  In'  tundamenta I  and  ode!  harmonic  s  arc 
mtmitesmal  at  tliis  time.  Since  only  the  tundamental  is  unstable-  tor  this  motor,  the  eye  i 
harmonics  dec  ay  with  time,  while  the  lmdament.il  begins  to  grow.  In  the  tune  inters  d 
from  alxiut  t  =  20  t < >  40,  lie  waveform  becomes  complex  as  it  transitions  tre.n  a  steep 
second  dominated  wave  to  an  almost  sinuso-dal  wave  at  the  tundac-i-nt  al  t  recjcieMc  s  . 
kecentlv,  tins  solution  was  continued  out  to  a  ninclimensioii.il  time  ot  Wto.  lhe  wave 
continues  to  decay  out  to  a  time  ot  aluntt  OH).  \t  this  time,  the  wave  amplitude  is  mils 
1.0,,  ot  the  mean  pressure  (compared  to  40,,  initially')  and  the  wave  is  essentially  a  pure 
fundamental  sine  wave.  \tter  t  —  100,  the  wave*  starts  to  grow,  and  is  still  growing  at 
t  =  100.  \s  the  wave  grows,  higher  harmonic  content  again  begins  to  .ippe.ir  as  a  result 
ot  energy  trarister  from  the-  tiindamental.  it  is  c-spe  ,ed  that  the  oscillation  will  continue 
to  grow  until  it  reac  lies  the  same  limiting  amplitndi-  as  the  other  test  cases  lor  this  motor 
and  propellant* 
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NONDIMENS  I0NA1,  TIME 


KONDIMENS I ONAL  TIJE 
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'■'!**  nii^m.il  jI  »|**i  l  i\t  * .  to  <  !rv»*l«  >p  .1  mr  f  ■  u  >< !  *  .  i  i » .  t  i  >  I « •  '»t  viKinj;  st rrp  -  t  roiit  r<;  r  »<  >  / 1 1  i  f  i  *  • ,  i  r 
i  »  M  .i  hi  i  it  \  1 1  »i  i| .[» -  i  .in  t  - 1  s.  1 1  okr  U  -s  s  t  .it  t  u  .i  1  rm  m  l  n  m  it  •  >  r  s  ,  I  i,i  s  i  in*n  .n  i  mi  ip  1 1  s|  iru. 

o  i  ’  uin  •  i  i*  .  t '  !,ii‘!'Iih!s  if  it  r  oilut  r  » *  r  r.  >  r  s ,  ,mil  t  hr  nm  lie  f  u  hi  met  h<  >o  hntl  \  hr  phs  x  u  h  I 

.  r.»i.i<-  i  i . .  i  m  •  m  •  l  vr  n  i  in  ‘i  1  In  i)»*  • 1  hi  t .  hru.  *  i  >  r  thr  prrxrnt  \,  r«  iblm  n,  ,i  mrt  h  >d  m.im-u  on  tin* 

<  t  j.  i  ih"t,i !  h  >i  i  .  ■  t  tin  ,i x  -  .*<•  in;  r-  *  1 1 ,  M\  ?>  r  n  j,  unu  \  r  1 1 1  u  i.i  I  <  omp  r»  ■  xxu  111  si  hrmr  x  \\  .is  found  t  ■  j 

'  “  -r.  >*  I  I>.  .!  t  t  l  hr  1  >1  h  r  r  xt  hrmrs  trslrd.  It  1ms  hrrn  shoAn  th.lt  tills  si  hrmr  is  t  hp.ihh- 

>  f  I*  ■  s.  r  i :  mi  i,,  i  s  hi  n  *«i  .is  .i  Si  i.i  r ,  >  him  nut  mints  w  it  1  h  nit  ■  rnr  r  it  in  ^  .irlitui.il  ; )  rr  -  .  >  t  ;  x  ist  - 

S.--H  k  is.  iliutiniis.  1  hi*  i  ■  i*  *t  *  u  >*  i  dors  not  rrl\  on  thr  use  .>1  am  .irtilu  i .  i  I  vistoxitv  .mo  is 

•  U  r'l*  it  r .  r  «*s<  •  r  v  !:ip  tin*  iii/.h  -  I  ro«|u«*tn  \  c  •  >i  it « 1 1  it  of  the  w.ivrtormx.  lhis  lomiuiihtion 
f  ■  •*  hn  i>  j  i  i*  •  t  .in  n  I  s< .  t  ff.it  t  hr  rr  t  Irt  t  ion  of  shot  k s  f  rom  i  *  nmo.i  r  u*  s  jruj  h«i  x  x;  i  m  1 1  dittuxur 

.ii":  - 1 1  >,  r  >  i \f  rff'ifs  «*v«*n  ntt**r  n.hn\  a. nr  iwlrs* 

I  )«•  i :  •  1 1  it  \  to  .UKilv/t*  sprttfiilU  tin*  1 1  MDput  **(i  results  A.is  <uldrd  tn  tin*  noitlinr.ir 
tost  iimitv  p  f- v,  r.i-ii.  Inis  i.ip.ilulitv  simplitifs  thi*  into  rpfft.it  ion  t»t  vomplrx  Ahvrtormx  .mo 

A  i  i  l  t.U  lilt  itr  i  4  >:  i  Ip  .1  f  I  soil  s  A  It  :t  Motor  dht.l  iiiuj  v\ltll  tilt*  rrxultx  oj  .ippruximhtr  nonlinrhr 

lot  Hh.. 

h-  .ntdil  iuii  to  thr  .impi  it  ndr,  Ahvrtorm,  .inti  lot  .it  ion  of  tnt*  prrsxurr  pulxr,  thr 
oft  nis  »t  thr  ismm  i.itotj  vt*|iw  it  v  (iist  iim.int  r  ,irr  jlxn  vrry  import,  int  in  drt  rrnnn  motor 
r  t  *  S|  m  i  f  i  si  ■  i'voi  aIm-m  v«|im  il\  <  oupl  in^  is  not  <  onsidrrcu. 

I ’u  solution  ot  tin*  mxt  ibil  it  \  problr  n  without  partu  Irs  ri'.u  find  .1  limitmj;  amplitude 
im  Jrprndrnt  ot  thr  initi.il  disturbance  .impi  itudr  .imp  A.ivrtorin.  I.irlirr 

<  .1 1<  ul  .it  ions  Aith  p.iftul<*s/  slmsxnl  limiting  .iniphtudo  to  hr  .1  st  ronj.;  tunition  ot 

tnr  miti.il  list  nr  i.uu  r.  Ihus,  it  .iium.irs  th.it  nonlmo.ir  [>,irtulc  (l.iuipinp  rtints  pins  ,t 
sp.,mtu  .1  it  r  i|»*  in  virt  *•  finniinn  thr  limit  it  1 1;  .implitml*’  ot  ,111  in  st  <  1 1  >  1 1  it  \  . 

I  hr  ns.-  ot  irlifit  l.l  I  v  1  s  <  o  s  1 1  \  insults  in  thr  nonph\sn.il  ,i  1 1  riui.it  ion  « >  t  thr  lu.’.h 

ti.iriiMih  s.  Inns,  t  i»r  use  ot  <i  mnn<  rit  .il  mrt  hoo  a  h  n  h  r  r  I  ir  s  on  ,m  .irtitui.il  v  im  os  it  \  to 


il.imp  [H»st-siioik  osi  illnt  ions  is  not  rr<  oip  nriniri!  tor  pmMr ms  fh.il  rrpuirr  thr 
»’xtrn»J  over  m.inv  a.u**  i\r  Irs. 


»i  ut  i<  >n  to 


I  Mil  J  '  )t  SYMBOLS 


c  r  is  s  sect  ji >s i  area 
nus  < in S' ,  sound  speed 

specific  lusit  ot  tL.i'.  at  constant  pressure 
siHS'itic  heat  of  sain)  pr. >()«''  ant 

par  t  ic  l<*  .»s  inter  a<  I  ion  Y.  >r  i  «■  [  >«  •  r  unit  volune,  t  ■  »r  tin 

gr.  nip 

'■las',  flux  t  mm  Bur  linn  s. it  ace 
number  if  particle  nr>Hips 
pressu  re 

particle  nas  heat  transfer  rate  per  unit  volune,  tor  t 
U  rixip 

t  unpe  r  atnre 

flame  temperature  of  the  pmpellant 
tenper ature  at  the  propellant  surface 
backwall  temperature  of  the  propellant 
t  ime 
veloc  i  t  v 

velocity  of  tin-  ounlmst  ion  pr.nlm  Is  as  thev  enter  the 
axial  component  of  us 
axial  distance 

total  particle  to  nas  weight  floss  ratio,  ,,  , 
particle  to  t;as  .veinht  flow  ratio  of  the  i(  1  particle  | 
nas  on V  isentmpie  exponent 
hens  itv 

mass  Uiroinn  rate,  per  unit  length,  per  unit  cross  set 


*  i1  n  part  ic  I  * 


he  |t'i  partic 


uia in  tiovs 


’.roup 


t  ional  area 


I  Mil  f  Of  SYVUSOI  S  (Continued) 
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Subscripts 

f  flame 

f! 

pj  i1*1  particle  group 

s  at  the  burning  surface 

t  at  the  nozzle  throat 
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